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RESONANCE TYPE INSTABILITIES 
IN THE GASEOUS DISKS 
OF THE FLAT GALAXIES 
II. The stability of solitary vortex sheet 

C.M.Bezborodov and V.V.Mustsevoy 

Linear stability analysis of the axisymmetric interface of velocity and density discontinuity 
in rotating gaseous disk has been performed numerically and analytically. Physical mechanisms 
leading to development of centrifugal and Kelvin-Helmholtz instability at the kink has been 
analysed in detail. In the incompressible limit it has been shown in the first time such areas in the 
parameter space that Kelvin-Helmholtz instability is stabilized by the density kink. This effect 
is caused by both specihcal angular momentum conservation and buoyancy. The possibility of 
application of obtained results to the stability analysis of the gaseous disks of the real flat galaxies 
is discussed. 

Introduction. 

The interest to a problem of stability of gaseous galactic disks is caused by the hydro- 
dynamic concept of an origin of spiral structure of galaxies proposed in 1972 by A.M. Frid¬ 
man (Fridman 1978). According to this concept the galactic spiral arms represent waves of 
density, arising and then reaching nonlinear stage as result of hydrodynamic (non-Jeans) 
instabilities development in a gaseous disk. 

The most probable applicants for a role of the spiral structure generator are thus 
Kelvin-Helmholtz and centrifugal instability. The development of these instabilities on a 
shift layer in rotating cylindrical and disk configurations of gas was already investigated 
many times (see Morozov et al. 1975, 1976a, 1976b, Morozov 1977, Landau 1944). In ac¬ 
cording to Morozov (1979) and Torgashin (1986) development of hydrodynamic instability 
in model with smooth VJ-like rotation curve (i.e. having significant velocity fall between 
solid-body rotating and plato parts) was considered, and the jump of angular speed was 
concurrent to jump of density. 

However, because of impossibility of analytical study the solution was found numer¬ 
ically, and from multiparametricity of a problem it was not possible to trace dependence 
of instability properties (growth rate) from parameters of model in details. Nevertheless 
we note that Morozov (1989) reported on instability weaking with growth of density con¬ 
trast between inner and outer disk areas and about minority of influence of sound speed 
contrast. In all the other quoted works the analytical asymptotic solutions were received 
in uniform density models 1 . 

We continue study of sharp speed jump stability, concurrent on radius with density 
jump, and promote it in discontinuous non-homentopic model. 

In the present part of the article we shall use the references to the formulae from a 
Part I, adding in this case symbol “1.” before number of formula. 

1 It is impossible to bypass work of Fridman (1989), where the deep analysis of gradient 
instabilities in a gaseous disk in the application to the hydrodynamic concept of spiral structure 
origin is spent. Unfortunately in the paragraph devoted to stability of velocity jump concurrent 
with density jump, faulty boundary condition (then corrected by Frirnan & Khoruzhii (1993)) was 
used, but investigation of growth rates of instability was not carried out. 
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1. Model and dispersion equation. 

The most simple and convenient, for the analysis model with density profile dis¬ 
continuity cannot be polytropic (Morozov & Mustsevoy 1989). It is connected that at 
p(r) — Jia 7 (r), where K = const, the density jump existence is equivalent to jump of 
pressure. The same takes place for three-dimensional quantities too. By virtue of a condi¬ 
tion of radial forces balance: 

Yl = l c K + c ll (i) 

r p dr dr 1 

where V v is gas rotation linear speed, T is gravitational potential, the pressure jump with 
necessity entails jump of gravitational potential (the potential of double layer). 

The data on distributions of galactic clouds velocities dispersion are practically 
absent; we shall use non-homentopic model with density jump but uniform pressure, 
that imposes relation on three-dimensional parameters of gas at both sides from jump: 
p ex / Pin — c sm/ c s e x (^ ie index “in” designates quantity from the inner side of jump and 
index “eaf with outer one). 

Considering told we research dynamics of disturbances of small amplitude on a back¬ 
ground of thin three-dimensional non-selfgravitating non-homentopic disk unlimited on 
radius being in external stationary axisymmetic potential and consisting of ideal noil- 
viscous gas. Following a standard linearization procedure we present all describing disk 
magnitudes as follows: 

/(r, <p, z, t ) = / 0 (r, z) + f(r, <p, z, t ), where |/| < f 0 . ( 2 ) 


The equilibrum axisymmetic stationary condition described by the radial balance equation 
( 1 ) is characterized by the following parameters distributions: 


Po(r, Z ) = Pin{z) + {Pex(z) - p in (z )) 9(r - R ), 


c s (r, z) = c sin (z ) 



1 0(r - R) 


Q(r, z) = ttin(z) + (^ex(z) - ^in(z)) 0{r - R ), 
Po(r, z) = c 2 s (r , z)p 0 (r, z )/y = const(z), 


( 3 ) 


where 9(x) is centered Heaviside function, 9(x) — 0 at x < 0, 9(x) — 1 at x > 0 and 
9(0) — 1/2, and 7 is adiabatic index. 

We do not use the Poisson’s equation determining gravitational potential because 
dealing with model problem we shall hereafter implicitly set necessary for us form of 
To (r, z) believing condition (1) for given distributions (3) and condition of hydrostatic 
balance along ^-coordinate holding. 

Because of uniformity of model parameters on azimuthal angle (p and time t the 
solution of gas dynamics linearized equations system is searched by decomposition on 
normal modes. Thus separate Fourier-harmonics of a kind / oc f(r, z) exp(im.(p — iut) 
evolution is further investigated. As Morozov (1989) show, for non-homentopic in a disk 
plane model (Po(r, z) /pg(r, z) ^ F(z )) it is described by the system of equations: 
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Here designations of a Part I are kept, i.e. p is disturbance of pressure, £ is perturbed radial 
Lagrange displacement, correlated to disturbance of a radial speed component by relation: 
d£/dt — —i(u> — mf2)£ = — iu)£ = fi r , cD is Doppler frequency, ft 2 = 20(20 + rdVt/dr ) 
is squared epicyclic frequency, c s is sound speed, the prime stands for differentiation on 
radius. 

The system (4)-(5) adequately describes dynamics of small disturbances in such ob¬ 
jects as: 

• a differentially rotating cylindrical configuration of non-dissipative non- 
selfgravitating ideal gas with the equation of state Po = po c s/7 (-Po and Po are homogeneous 
along cylinder forming equilibrum pressure and density respectively, c s is adiabatic sound 
speed, 7 is adiabatic index); the description is adequate for disturbances with v z = 0 ; 

• a thin non-dissipative non-selfgravitating disk of polytropic gas, being in stationary 
external (star) gravitational potential of a kind T(r, z) = <F(r) + x( z ) ( see Appendix I in 
the work of Gorkavy & Fridman (1994)); pinch type disturbances on ^-coordinate are in 
this case considered only, v r and v^ with necessity designates their average on z values, 
and Po and Po are surface pressure and density accordingly; the two-dimensional adiabatic 
index y s is associated to a volumetric by relation: y s = 3— 4 /( 7 +1) (Churilov & Shukhman 
1981). 

Note that the realization of the last case is necessary for construction of step- 
homogeneous two-dimension model of thin disk because for potential of a more general 
form T = T(r, z) either rotation velocity averaging V v = V^(r, z) on ^-coordinate will not 
lead to discontinuous distribution of V^(r) or there will be the radial inhomogeneities of 
pressure, which are necessary for fulfilment of equilibrum forces balance. A really specified 
form of potential means that the disk is quite thin for applicability of potential decom¬ 
position in a series on small parameter z/r with neglecting of all terms more senior than 
second-order. 

As in each of areas divided by jump all equilibrum parameters are homogeneous, the 
system (4)-(5) becomes identical to system (1.1)-(1.2) and is reduced to modified Bessel 
equation on p, the general solution of which is a linear combination of modified Bessel 
functions I m {kir ) and K m (kir ) where a designation is introduced 


in\ex 


^^in;e x m ^in;ex) 

c 2 . 
sin;ex 


( 6 ) 


Note that completely analogous equation is also valid for quantities p(r), rj(r) — p(r)/po, 
V (r) = c 2 s r]{r). Assuming hereafter for uniqueness of the solution Reki n - ex > 0, from a 
natural condition of its limitation in a disk centre and on infinity follows: 


f (/CjfjT"), T P, 

\BK m (k ex r), r > R. 


( 7 ) 


Boundary conditions for the solutions on jump we determine by integrating the equa¬ 
tion (4), (5) on radius from R — £ up to R + £ with account of peculiarity of chosen model 
- dPo/dr — 0 and then tending £ —> 0 (see Fridman & Horuzhii 1993 for details): 


£ (P + 0) — £ (R — 0) = 0, 


( 8 ) 
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( 9 ) 


P (R +o) - P (r - o) = Pi ~ + ((r) ml. ~ ^L). 

Coupling of the solutions at r — R results system of the equations on factors A and B. 
The condition of simultaneity of this system represents the dispersion equation allowing 
to determine frequency eigenvalues: 


Pin Jpk‘i n R 2 ( Xe X 

1 k 2 ex R 2 -^(l + Q)(l-q 2 )a. 


( 10 ) 


Here the following designations are made: M = R il ir , /c srn is Mach number, q — 
ttex/ftin, Q = Pin/pex, P = c sex fc ain , and /3j are defined by (1.8) and (1.9) in which 
the index i still takes values “in” and “ex”, and the top index is omitted because of 
coincidence of radii of speed and density jumps. Besides let us introduce dimensionless 
frequency x — co/fli n . The existence of positive Imuo (growth rate) means instability. 

Note the original dispersion equation was satisfied identically for neutral values 
(Imuo = 0): to — (m±2)O in , uo — (m±2)Tl ex , that is for gyroscopic modes of fluctuations 
of a homogeneous solid-body rotating disk. It is natural result, as the disturbances with 
such frequencies propagate exclusively in an azimuthal direction (/c in = 0, or k ex — 0) 
and are not “sensitive” to radial inhomogeneities. Accordingly the given modes do not 
influence system stability and deriving (10) we supposed that the frequency does not take 
such values. 

2. The stability analysis. 

Let us obtain the approximated solutions of the dispersion equation (10) for limiting 
cases of large and small compressibility of medium. 

a) The case of large compressibility (.M S> 1) 

By applying representation of arguments of Bessel functions through dimensionless 
parameters note the module of argument is dictated first of all by Mach number M. 
Assuming therefore \ki n R.\ 3> to, \k ex R\ 3> to, from (10) we find: 


x - n + q + q>1 ) + ^ M (! + ( 5)( 1 ~ q2 )| • ( n ) 

See |x| ~ M, \k ex R\ ~ M 2 , ~ M 2 , therefore the assumptions made earlier are 

justified backdating. 

For a case of constant density from (11) we get result obtained by Morozov (1977): 

uo ~ [ m (i + q) + iM( 1 — q 2 )] . (12) 

2 


Taking into account that expression // = \/Q takes place in model without self-gravitation 
let us write down growth rate (11) as: 

( 13 > 

Remarkable peculiarity of obtained result is extremely weak dependence of growth rate on 
magnitude of density jump at pi n f> (:X . and on the contrary proportionality to 1 1 at 
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Pin -C p e x- Applying similar reasons we determine that for wave pattern angular rotation 
phase speed fi p = Recu/m the following estimations fi p ~ at pi n p ex and fi p ~ Vt ex 
at pi n < pe* are valid. Hence the wave pattern rotation speed comes nearer to rotation 
speed of more dense medium at significant density differences. 

Note that in absence of speed shift (q = 1) neutral solution u — mVti n respecting 
to 11011 -growing disturbances “frozen” in gas follows from (11). Thus density kink alone 
cannot cause instability. 

To shed some light on physics of instability it is convenient to pass from (11) to the 
dimensional form: 


u ~ 


A7 (pin^s in^in + Pex^sex ^ea:) t 2 (P in + p ex )mi. - np 
CsinPin T C sex p ex 


(14) 


From (14) follows that the growth rate of considered instability in the main order is equal to 
the ratio of difference of density of centrifugal force, acting on a liquid particle fluctuating 
in small vicinity of jump, to the sum of wave resistance (characteristic impedances) of 
media on both sides from jump. Thus this instability must be called centrifugal one, just 
same as in case of disk of uniform density, which Morozov (1977) had consider 2 . 

b) An incompressible limit (M <C 1). 

Taking into account \k vn - ex R\ oc M we make following assumption: \ki n R\ <C 1, 
\k ex R\ <C 1. Applied asymptotic representation of modified Bessel functions at small 
values of argument (Handbook of Mathematical Functions... 1964), we find the solution 
(10) in the dimensionless and dimensional forms: 


x ~ 


1 

1 + Q 


| m(q + Q) + {q-Q) ± 


± 


(< q~Q ) 2 


im(l-Q) 2 (l-g 2 ) 


m 2 Q( 1 



(15) 


tO — % ^{yiinPin T H exPex ) T ( ^exPex H inpin ) T 

Pin T Pex t 

^{^SlexPex ^in Pin ) — Tn(p ex Pin) (^*n ^ea;) 



The negativity of expression under square root sign in (15) and (16) means instability. 

In a case <5 = 1 (uniform on density at z — const disk) the result received by Morozov 
(1977) follows from (15) : 

u ~ | m(q + 1) + (q — 1) + i\l — q\\Jm 2 — l| . (17) 

2 Term “centrifugal” with reference to considered instability was introduced later (Morozov et 
al. 1984, 1985a, 1985b). 
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In absence of velocity difference on jump (q = 1) neutral solutions follow from (15) as well 
as in the case of large compressibility, but in this case they are the pair: u — rnii tri and 

u ~ Cl in + 21 + ^ ■ (18) 

At significant density contrast (18) governs gyroscopic modes: u ~ (to + 2 )fi in at Q -C 1 
and u; ~ (m — 2)Q m at Q 1. As one can see at M <C 1 the density jump alone cannot 
cause instability too. 

Let us analyse a physical nature of terms in (15) and (16) which stabilize and desta¬ 
bilize the system. 

For the disturbances with shortest azimuthal wavelength (m 3> 1) 


U ~ TO- 


^ex Pex A inPi 
Pin + P ex 


%Vn 4 / pinpex 


iS^in ^ex)^ 
{Pin "1” Pex ) 2 


(19) 


follows from (16). 

Introducing azimuthal wave number on jump k v = m/R and linear speeds of rotation 
inside and outside in immediate jump proximity Vj„ = /?Q m and V ex = accordingly, 

it is easy to bring (19) to the form: 


, Pex^ex T Pin^in , - 7 V PinPex i 

W ~ -—--- \V in -V ex \. (20) 

Pin T Pex Pin T Pex 

The last expression is identical to classical expression describing Kelvin-Helmholtz insta¬ 
bility (KHI) developing on vortex sheet in plainparallel current of an incompressible fluid 
in absence of gravitation force (g = 0): 


, Pi hi + P2V2 
u = k 11 -± 

Pi + P 2 



P 2 ~ Pi 
Pi + P2 


, 2 Pip 2 (Vi-V2) 2 
11 (P 1 + P 2) 2 


( 21 ) 


Here the index 1 marks magnitudes concerning to the top layer, 2 to the bottom one respec¬ 
tively, k\\ is wave number of disturbances along vortex sheet 3 . Thus the nature of the last 
term under square root sign in (16) always destabilizing jump surface is perfectly trans¬ 
parent: it is caused by Bernoulli effect. This result is quite expected as for disturbances 
with short on cp wavelength effects of curvature and rotation are the least essential. 

Second term under square root sign in (16) can be rewritten as follows: 

_7 (Pex ~ Pin ) R^in ~ R^ex (OOl 

Pex + Pin ) 2 2 1 J 

This term is directly proportional to centrifugal acceleration difference on jump, renders 
stabilizing influence at fl ex > 0* n and destabilizing at fl ex < Qj n . So the effects of 
centrifugal stabilization and destabilization (see, for example, Nezlin & Snezhkin 1990) 

3 Inasmuch as our model is quasi-two-dimensional (v z = 0), equation (21) is written out in 
the assumption of concurrence of directions Vi, V 2 and kn. 
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take place, which are well known for a compressed fluid. However, for an incompressible 
liquid they are shown for the first time. It is also interesting that these effects act only in 
a case pi n ^ p ex . It has resulted to such fact that they were not found by Morozov (1977). 
The latter, on our sight, is a direct consequence of incompessibility and is explained as 
follows. 

The incompressible fluid is characterized by significant (infinite as an ideal) elasticity, 
because of that centrifugal effects cannot result instability in the same way, as in compress¬ 
ible fluid, where liquid particle releasing from vortex sheet by non-compensated centrifugal 
or gravitational force is displaced on radius “tightening” particles being ahead it. In given 
case the radial rearrangements of liquid particles can come true only by their mutual re¬ 
placement against each other, as it occurs at development of Rayleigh-Taylor instability, 
when particle of heavy fluid “sinks” on a place of an emerging particle of light fluid. Not 
casually therefore similarity of a considered combination (22) with first term under square 
root sign in (21), just governing the development of Rayleigh-Taylor instability on liquids 
interface in plainparallel flow. This similarity becomes especially transparent in the case of 
significant density contrast: assuming in (21) pi S> P 2 we get k^g and in the case p\ -C p-j 
this term gives —k\\g. From (22) as at pi n S> p ex and at pi n <C p ex follows 


-fee 




ml 


(23) 


Thus in both cases we have longitudinal wave number and mass force acceleration on vortex 
sheet production. Told leads to unexpected conclusion that for considered axisymmetic 
vortex sheet a combination of parameters 


Pin Pex (24) 

Pin + P ex 2 

is equivalent g in a plainparallel case. 

Generalizing told we conclude that (22) describes buoyancy effects in considered model 
and it is caused by buoyancy forces disbalance. 

At last the first term under square root sign in (16) is easily represented as 


|PexR fd ex pi n R 
PexR T PinR 


(25) 


which is possible to interpete as the ratio of a half-difference of density of an angular 
momentum on both jump sides to specific (per volume unit) momentum of inertia of gas 
on jump. Tracing genesis of this term in any case rendering stabilizing influence upon 
jump it is possible to note that it is caused by action of Coriolis’ forces and reflects the 
tendency to preserve an angular momentum of system. 

Rather interesting result is appearence of stability area on parameter q in a vicinity 
of q = 1 at Q ^ 1. This area wider as stronger Q differs from unit and less as the less is 
mode number to. The borders of stability constructed on the asymptotic formula (15) for 
to = 2, 3, 4 are shown on Fig. 2.1. 

On Fig. 2.2-2.5 results of the numerical solution of the dispersion equation (11) by 
iterative Newton-Rafson method are drawn. From them follows that the approximated 
formula (15) will be good coordinated with dispersion curves in an incompressible limit. 
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The same good consent of results with (11) takes place in a inverse case of essentially 
supersonic speed difference on jump, if both M 2 ( 1 — q 2 ) S> to and M 2 (l — q 2 ) / g 3> to 
simultaneously are valid. 

On our sight Fig.2.5 is rather interesting, on which stabilization of centrifugal mode 
with growth q is shown. Though it is occurs formally at M > 1 the point of stabilization 
and branching on two neutral modes corresponds to subsonic speed difference on jump 
M(1 — q)/n — 0.3 < 1, therefore the mechanism of stabilization is similar to described 
above. 

Thus comparative analysis of results of numerical and asymptotic analytical researches 
of the equation (10) allows to conclude that all conclusions made in this section are valid. 

3. Conclusions. 

Let us summarize the basic conclusions concerning to stability of non-homentopic 
angular speed and density discontinuity profiles of finite thickness (£ -C A^, A p , where 
An, A p are characteristic radial scales of angular speed and density jumps) in model of 
uniform on pressure at z — const gaseous disk. 

1. At essentially supersonic speed difference on jump centrifugal instability is raised, 
which growth rate is equal to the ratio of difference of density of centrifugal force acting 
on liquid particles in vicinity of jump, to the sum of media wave resistances at both sides 
of jump. Thus the angular rotation speed of a spiral pattern in case of significant density 
drop is close to more dense medium rotation speed. 

2. At small compressibility of both media the Kelvin-Helmholtz instability develops. 
It is caused by Bernoulli effect and at presence of density difference also by buoyancy effects, 
for which the role of external mass force density is played by the following combination of 
parameters: 

_ Pin ~ Pex ~ ex 

^ Pin + Pex 2 

Depending on sign of g e ff, an additional centrifugal destabilization or stabilization of 
KHI takes place. Besides, at comparison with KHI on plainparallel vortex sheet there 
are the effects caused by angular momentum preservation. This effects exert in any case 
a stabilizing influence, the more essential, the larger is disturbance scale (and the less 
number m) and the more considerable density drop. 

3. Density drop alone in absence of speed difference is not capable to result instability 
of considered model both in supersonic and subsonic case. 

It is necessary to note that as the KHI does not develop for disturbances with wave¬ 
length comparable to or exceeding characteristic scale Aq (that is modes with m > 1 
will appear neutral in view of finite jump thickness) and the large-scale disturbances are 
stabilized by virtue of conclusion 2, for smooth subsonic speed jumps there is the area 
of stability on parameters q, Q. The latter is necessary to take into account at stability 
analysis of real objects . 

The stabilization of Kelvin-Helmholtz instability for long wavelength (with small to) 
disturbances of subsonic non-homentopic rotation speed and density jump in view of sup¬ 
pression of this instability for short-wave disturbances because of inevitably present finite 
radial “smearing” of jump means an opportunity of rather long-lasting existence of such 
jump. Thus the rotation speed on radius of jump can be even supersonic. Such jumps 
are observed in central areas of some flat, galaxies (Afanesiev et al. 1988a, 1988b, 1988c, 



1991, 1992) and can play a role of an internal reflecting surface for waves amplifying on 
a corotation resonance (Fridman et al. 1994). Our analysis allows to assert that in case 
of significant density contrast on jump the life-time of such surface will be sufficient for 
development of resonant instability of the specified waveguide layer. Otherwise the jump 
will be smeared out because of “faster” surface mode development earlier than resonant 
disturbances will grow to nonlinear amplitudes. 
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Fig.2.1. Border of marginal stability of azimuthal modes m = 2. 3. 4 (accordingly, 
continuous line, long dash, short dash) according to the asymptotic formula (15). 
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Fig.2.3. Dependence of dimensionless frequency Reuo/fli n (a), (c) and dimensionless 
growth rate Imu)/tti n (b), (d) from Mach number on the data of the numerical solution 
(10) (continuous lines) and asymptotics (11) (dashed lines) for various sets of parameters 
m, q, Q. 
















